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Inclusive jet cross sections in Z/^* events, with Z/^* decaying into an electron-positron pair, 
are measured as a function of jet transverse momentum and jet multiplicity in pp collisions at 
= 1-96 TeV with the upgraded Collider Detector at Fermilab in Run II, based on an integrated 
luminosity of 1.7 fb~^. The measurements cover the rapidity region Ij/-""'! < 2.1 and the transverse 
momentum range pl^ > 30 GeV/c. Next-to-leading order perturbative QCD predictions are in good 
agreement with the measured cross sections. 

PACS numbers: 12.38.Aw, 13.85.-t, 13.87.-a 
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The measurement of the inclusive production of col- 
limated jets of hadrons in association with a Z j^* bo- 
son in "pp collisions provides a stringent test of pertur- 
bative quantum chromodynamics (pQCD) The un- 
derstanding of Z/7*-|-jets final states is a crucial part 
of the physics program at the Tevatron since they con- 
stitute important irreducible backgrounds in searches for 
new physics. Previous results [2] from Run I at the Teva- 
tron have been compared with leading-order (LO) plus 
parton shower Monte Carlo predictions affected by large 
scale uncertainties. This Letter reports new and more 
precise measurements of the inclusive jet cross sections 
in Z/^*{-^ e+e") production using 1.7 fb"""^ of data col- 
lected by the CDF experiment in Run II. Inclusive jet 
differential cross sections as a function of jet transverse 
momentum pi^' and total cross sections as a func- 
tion of jet multiplicity A'jet are measured. The data are 
compared to next-to-leading order (NLO) pQCD predic- 
tions [3] including non-perturbative contributions. 

The CDF II detector is described in detail else- 
where 0. The detector has a charged particle track- 
ing system immersed in a 1.4 T magnetic field aligned 
coaxially with the beam line that provides tracking cov- 
erage in the pseudorapidity Q range \ri\ < 2. Segmented 
sampling calorimeters, arranged in a projective tower ge- 
ometry, surround the tracking system and measure the 
energy of interacting particles for \ri\ < 3.6. The central 
electromagnetic and hadronic calorimeters 6] cover the 
region |?7| < 1, while the end- wall hadronic calorimeter Q 
provides forward coverage out to |ry| < 1.3. Forward elec- 
tromagnetic and hadronic calorimeters Q cover the re- 
gions 1.1 < I77I < 3.6 and 1.3 < |?7| < 3.6, respectively. 
The calorimeters are instrumented with finely segmented 
detectors jgl, Q to measure the shower profile at a longi- 
tudinal depth close to the location of a typical electro- 
magnetic shower maximum. Cherenkov counters in the 
region 3.7 < I77I < 4.7 measure the number of inelastic pp 
collisions to compute the luminosity [l^] . 

Monte Carlo event samples are used to determine de- 
tector acceptance and reconstruction efficiency, estimate 
background contributions, and unfold the measurements 
back to the hadron level [ll|. Samples of simulated 
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inclusive Z/j*{^ e+e^)-|-jets events have been gener- 
ated using the pythia 6.216 12] Monte Carlo generator. 
CTEQ5L [3] parton distribution functions (PDFs) are 
used for the proton and antiproton. The pythia sam- 
ples have been created using a special tuned set of pa- 
rameters, denoted as pythia-tune a 14], that includes 
enhanced contributions from initial-state gluon radiation 
and secondary parton interactions between proton and 
antiproton beam remnants and provides an accurate de- 
scription of the measured jet shapes and energy flows 



in Z/7*(— » e"'"e~)-|-jets final states [15|. Monte Carlo 
samples for background processes are generated using 
pythia-tune a. The samples are passed through a full 
CDF detector simulation (based on GEANt3 [16|] where 
the GFLASH [13] package is used to simulate the energy 
deposition in the calorimeters) and reconstructed and an- 
alyzed with the same analysis chain as for the data. 

Events are collected using a three-level trigger sys- 
tem [III- At the first-level trigger, events are required 
to have a central electromagnetic calorimeter cluster 
(l?7l < 1) with Et [3] above 8 GeV and an associated 
track withp^^'^'^ above 8 GeV/c. Similarly, at the second- 
level (third-level) trigger a central electromagnetic clus- 
ter with Et > 16 GeV {Et > 18 GeV) and an asso- 
ciated track with p*^'"^^ > 8 GeV/c {pt^'"=^ > 9 GeV/c) 
are required. The events are then required to have two 
electrons [3l with E^ > 25 GeV and a reconstructed 
invariant mass in the range 66 < M^e < 116 GeV/c^ 
around the Z boson mass. The electron candidates are 
reconstructed using criteria described in In this 

study, one electron is required to be central (jTy^J < 1) 
and fulfill tight selection cuts, while the second elec- 
tron is required to pass a looser selection and to be 
either central (CC final-state configuration) or forward 
with 1.2 < jry^l < 2.8 (CF final-state configuration). 
The trigger efficiencies for CC and CF configurations are 
99.96 ± 0.01% and 97.9 ± 0.3%, respectively. The events 
are selected to have a reconstructed primary vertex with 
z-position within 60 cm around the nominal interaction 
point, and at least one jet with corrected transverse mo- 
mentum p^x'cor ^ GeV/c (see below), rapidity [3] in 
the range ]?/^°j] < 2.1, and ARe-^jet > 0.7, where ARe-jet 
denotes the distance {y — 4> space) between the jet and 
each of the two electrons in the final state. The final 
sample contains 6203, 650, 57, and 2 events with at least 
one, two, three, and four jets, respectively. 

Jets are reconstructed in data and Monte Carlo simu- 
lated events from the energy deposits in the calorimeter 
towers with transverse momenta 2lJ above 0.1 GeV/c. 
The towers associated with the reconstructed electrons in 
the final state are excluded from the jet search. Jets are 
searched for using the midpoint algorithm 22] with cone 
radius R = 0.7 and a merging/splitting fraction of 0.75, 
starting from seed towers with transverse momenta above 
1 GeV/c. The same algorithm is applied to the final state 
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particles in the Monte Carlo generated events, excluding 
Z/7* decay products, to define jets at the hadron level. 

The rapidity and azimuthal angle of the jets, and 
"^^caii wG^l reconstructed in the calorimeter with a res- 
olution better than 0.05 units in both y and (p. The 
measured jet transverse momentum p^x'cai systematically 
underestimates that of the hadron-level jet. For p^x*cai 
values about 30 GeV/c, the jet transverse momentum 
is underestimated by about 30%. The systematic shift 
decreases with increasing pl^^^^i down to about 11% for 

Pr^cai ^ GeV/c. This is mainly attributed to the 
presence of inactive material and the non-compensating 
nature of the calorimeter [2^. An average correction, 
as a function of Px*cai ^^^1 y^^^i, is applied to the mea- 
sured p^x*cai account for these effects ji^]. The mea- 
sured P^x*cai ^^^'^ includes contributions from multiple pp 
interactions per crossing at high instantaneous luminos- 
ity. Multiple interactions are identified via the presence 
of additional primary vertices reconstructed from charged 
particles. For each jet, p^x'cai corrected for this effect 
by removing a certain amount of transverse momentum, 
i5™ — 1.06 ± 0.32 GeV/c, for each additional primary 
vertex in the event, as determined from data [24]. 

The main backgrounds to the Z/j*{~^ e^e^)-|-jets 
sample arise from inclusive-jets and M^-|-jets events, and 
are estimated from the data. First, an inclusive jet data 
sample is employed to estimate the probability for a 
jet to pass a given electron selection. The probabilities 
are parametrized as a function of ?5x*cai ^'^^ typi- 
cally around 0.02 and 0.005 for central and forward loose 
electrons, respectively. Second, a sample of events in 
data with exactly one reconstructed tight central elec- 
tron is selected. For each jet in the event, the of 
a fake electron is determined, and the invariant mass 
of the tight-central electron and jet final state is then 
computed. Event-by-event, all electron-jet combinations 
that fulfill the cuts and with an invariant mass within 
66 < Ale-jet < 116 GeV/c^ are considered in the back- 
ground calculation, where each combination is weighted 
by the corresponding value and divided by the num- 
ber of accepted electron-jet combinations in the event. 
The total inclusive-jets and iy-|-jets background is then 
computed in each measured distribution. Other back- 
ground contributions from tt, Z/j*{~^ e+e^) -t- 7, WW, 
WZ, ZZ, and Z/'j*{^ r+r^)+jets final states are esti- 
mated using Monte Carlo samples. The total background 
in inclusive Z/j*{-^ e+e~)-|-jets production is about 12% 
for A'jot > 1, and increases up to about 17% for Nj^t > 3. 
Good agreement is observed in the total number of events 
between the data and the Z/7*(— > e"'"e~)-|-jets signal plus 
background predictions. A test, where only statistical 
uncertainties are considered, gives probabilities that 
vary between 80% and 25% as iVjot increases. 

Raw inclusive jet differential cross sections as a 



function of Px'cor defined as da/dpl^^^^^ ~ Z ^ 
(-^jet'^/APT*cor)' whcre N^°^ denotes the total number of 
jets in a given p^^cor ^^^^ ^^'r^cor ^he size of the bin, 
and C is the luminosity. N^°^ is corrected bin-by-bin 
for background contributions and trigger inefficiencies. 
The measured cross sections are then corrected for accep- 
tance and smearing effects back to the hadron level using 
PYTHIA-TUNE A Monte Carlo event samples, and a bin- 
by-bin unfolding procedure that also accounts for the effi- 
ciency of the Z/'y*{-^ e+e^) selection criteria. The final 
results refer to hadron level jets with p!^^ > 30 GeV/c and 
\y'°^\ < 2.1, in a limited and well-defined kinematic range 
for the Z/7* decay products: E^ > 25 GeV, \r]'''^\ < 1.0, 
Irj^^l < 1.0 or 1.2 < Irf^j < 2.8, 66 < Mee < 116 GeV/c^, 
and ARe-jet > 0.7. In order to avoid any bias on the cor- 
rection factors due to the particular PDF set used, which 
translates into slightly different simulated p^x*cai distribu- 
tions, the PYTHIA-TUNE A Moutc Carlo event sample is 
re- weighted until it accurately follows the measured p^x'cai 
spectra. The unfolding factors Uip!^^^) = -^/ — — 

are computed separately for the different measurements 
and vary between 2.0 at low pi^^ and 2.3 at high p!^^ . 

A detailed__study of the systematic uncertainties was 
carried out 115]. A ±1.5% uncertainty on the trigger ef- 
ficiency translates into ±1.5% and ±0.06% uncertainties 
on the cross sections for CF and CC configurations, re- 
spectively. The uncertainty on the p!^^ dependence of the 
electron identification efficiency introduces a ±5% uncer- 
tainty on both CC and CF results. The measured jet en- 
ergies are varied by ±2% at low pi^^^ed to ±2.7% at high 
Pr^cai to account for the uncertainties on the absolute 
energy scale in the calorimeter this introduces un- 
certainties on the final measurements which vary between 
±5% at low pj^* and ±12% at high pl^\ The yj^* depen- 
dence of the average correction applied to P^^^^i intro- 
duces a ±2% uncertainty on the measured cross sections, 
approximately independent of p!^*". The uncertainty on 
(5™ has a negligible effect on the measured cross sec- 
tions. The uncertainty on the Px*;,^! dependence of /j"* 
introduces a ±15% uncertainty on the inclusive-jets and 
M^+jets background estimation, that translates into a 
less than 2% uncertainty on the measured cross sections. 
A conservative ±30% uncertainty on the normalization 
of the rest of the background contributions, as extracted 
from Monte Carlo samples, introduces a less than 1% ef- 
fect on the final results. If the unfolding procedure is 
carried out using unweighted pythia-tune A, the effect 
on the measured cross sections is less than 1%. Posi- 
tive and negative deviations with respect to the nominal 
cross section values are added separately in quadrature 
to define the total systematic uncertainty. The final re- 
sults are obtained from the combination of CC and CF 
measurements. Finally, a 5.8% uncertainty on the total 
luminosity is included in the measured cross sections. 
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FIG. 1: (a) Measured inclusive jet differential cross section 
as a function of pi^* (black dots) in Z/'y*{^ e^e~)+jets with 
A'jot > 1,2 compared to NLO pQCD predictions (open cir- 
cles). For clarity, the measurement for A'jot > 1 is scaled up 
by (x20). The shaded bands show the total systematic un- 
certainty, except for the 5.8% luminosity uncertainty, (b and 
c) Data/theory ratio as a function of pi^* for A'jot > 1 and 
A'jot > 2, respectively. The dashed and dotted lines indicate 
the PDF uncertainty and the variation with n of the NLO 
pQCD predictions, respectively. 



Figure [TJa) shows the measured inclusive jet dif- 
ferential cross sections as a function of pl^ 



e )-l-jets production, with iVjot > 1 



m 
and 
The 



-^jot > 2, compared to NLO pQCD predictions, 
data are reported in Table L The cross sections decrease 
by more than three orders of magnitude as pl^^ increases 
from 30 GeV/c up to about 300 GeV/c. The NLO pQCD 
predictions are computed using the MCFM program [4] 
with CTEQ6.1M PDFs 



PYTHIA-TUNE A Monte Carlo samples, as the ratio be- 
tween the nominal p!^^ distribution and the one obtained 
by turning off both the interactions between proton and 
antiproton remnants and the string fragmentation in the 
Monte Carlo samples. The correction decreases as p!^^ in- 
creases from about 1.2 (1.26) at pl^^ of 30 GeV/c to 1.02 
(1.01) for pi;;* > 200 GeV/c for TVjet > 1 (A^jet > 2), and 
is dominated by the underlying event contribution. In or- 
der to estimate the uncertainty on Chad, pythia samples 
are generated with a different set of parameters, denoted 
as TUNE DW [2^, that governs the underlying event ac- 
tivity and also describes the Z/j*{-^ e+e^)+jcts final 
states. The uncertainty on Chad is about 10% (17%) at 
low pi^ and goes down to 1% at high pi^ for Njet > 1 
(-^jct > 2). The ratios between data and theory as a func- 
tion of p^^;* are shown in Fig. [Hb,c). Good agreement is 
observed between the measured cross sections and the 
nominal theoretical predictions. A test, where the 
sources of systematic uncertainty on the data are con- 
sidered independent but fully correlated across p'r^^ bins, 
and the uncertainty on Chad is also included, gives a 
probability of 99% (22%) for Nj^t > 1 (^jct > 2). 
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with the renormalization 
and factorization scales set to /i^ = A/| + p'^{Z), and 
using a midpoint algorithm with R ~ 0.7 and i?scp — 
1.3 [2a] to reconstruct jets at the parton level. Values 
for i?sGp between 1.0 and 2.0 change the theoretical pre- 
diction by less than 2%. A variation of /z by a factor 
two (half) reduces (increases) the theoretical predictions 
by 10% to 15%. The uncertainties on the NLO pQCD 
predictions due to the PDFs were computed using the 

Hessian method . They vary from ±4% at low j^r^^ to xhe measured event 



FIG. 2: (a) Measured total cross section for inclusive jet pro- 
duction in Z/7*(— > e^e~) events as a function of A'jot com- 
pared to LO and NLO pQCD predictions. The shaded bands 
show the total systematic uncertainty, except for the 5.8% 
luminosity uncertainty, (b) Ratio of data and NLO to LO 
pQCD predictions versus A'jot. The dashed and dotted lines 
indicate the PDF uncertainty and the variation with fi of the 
NLO pQCD predictions, respectively. 



Figure [5] shows the total cross 
for Z/j*{-^ e"'"e~)-|-jets events up 



sections, 



to m 



jot J 



JCt 



> 



±10% at highp^T • 

The theoretical predictions include parton-to-hadron 
correction factors, Chad(A^jot,?3x*)' ^^^^ approximately 
account for non-perturbative contributions from the un- 
derlying event and fragmentation into hadrons (see Ta- 
ble I). In each measurement Chad is estimated using the 



cross sections are: ui = 
fb, (72 = 695 ± 
fb, and 0-3 = 60 ± 



7003 ± 146(stat.)+4?^(syst.) ± 406(lum 
37(stat.)t^^(syst.) ± 40(lum.) 
ll(stat.)tg(syst.) ± 3.5(lum.) fb, for Nj^t > 1, A^jot > 2, 
and iVjot > 3, respectively. The data are compared to 
LO and NLO pQCD predictions. The parton-to-hadron 
non-perturbative corrections vary between 1.1 and 1.4 as 
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Tel 

[GeV/c] 



± (Stat.) ± (syst.) ± (lum.) 

dpi. 



[fb/(GeV/c)] 



Chad ± (stat.) ± (syst.) 
parton — > hadron 



-)+jetS (JVjet > 1) 



30 - 


35 


413.3 ± 13.3+3'i''3 ± 24.0 


1.209 ±0.010 ±0.134 


35 - 


41 


263.3 ± 9.4l^^;| ± 15.3 


1.146 ±0.010 ±0.096 


41 - 


47 


178.3 ± 7.5tll ° ± 10.3 


1.114 ±0.011 ±0.077 


47- 


54 


128.5 ± 5.9tll ± 7.5 


1.097 ±0.012 ±0.066 


54- 


62 


80.5 ±4.3t^;^ ±4.7 


1.086 ±0.013 ±0.059 


62 - 


72 


52.5 ± 3.2t4 3 ± 3.0 


1.078 ±0.013 ±0.053 


72 - 


83 


34.2 ± 2.4l2;8 =•= 2-0 


1.072 ± 0.015 ± 0.049 


83 - 


110 


16.0 ±l.lli;3± 0.9 


1.063 ± 0.012 ± 0.043 


110 - 


146 


4.9 ± 0.5j:°-j ± 0.3 


1.051 ± 0.012 ± 0.035 


146 - 


195 


l.l±0.2t°;J ±0.06 


1.040 ± 0.008 ± 0.027 


195 - 


400 


0.08 ± 0.03+°;°^ ± 0.005 


1.021 ± 0.005 ± 0.013 



-)+jets (iVjet > 2) 



30 


- 38 


52.9 ± 3.5+4 g ±3.1 


1.262 ±0.022 ±0.217 


38 


- 47 


37.0 ±2.81:2.8 ±2.1 


1.207 ±0.024 ±0.169 


47 


- 59 


21.2 ±1.8li;|± 1.2 


1.164 ±0.025 ±0.130 


59 


- 79 


10.5 ±1.0t?;g± 0.6 


1.123 ±0.024 ±0.093 


79 


- 109 


5.7±0.6t°-j ±0.3 


1.087 ±0.026 ±0.062 


109 


- 179 


0.88±0.15lo;'j'g ±0.05 


1.052 ± 0.020 ± 0.030 


179 


- 300 


0.15 ± 0.04 1° °2 ±0.009 


1.026 ± 0.010 ± 0.008 



TABLE I: Measured inclusive jet differential cross section 
in Z/'y*{—> e'''e~)+jets production as a function of pi^* with 
A'jet > 1 and iVjet > 2. The systematic uncertainties are fully 
correlated across pi^* bins. The parton-to-hadron correction 
factors Chudip'T^ , Njet) applied to the pQCD predictions. 



Njet increases. For Njet > 1 and A'jet > 2, the LO pQCD 
predictions underestimate the measured cross sections by 
a factor about 1.4, which corresponds to probabihties 
of 0.07% and 2.7%, respectively. Good agreement is ob- 
served between data and NLO pQCD predictions, with 
probabilities better than 83%. For iVjet > 3, where no 
NLO pQCD prediction is available, the measured cross 
section indicates that the data can be described by a LO- 
to-NLO theoretical factor independent of A'jet. 

In summary, we report new results on inclusive jet 
production in Z/j*{-^ e+e") events in pp collisions at 
y/s = 1.96 TeV for jets with pl^* > 30 GeV/c and 
lyj''*! < 2.1, based on 1.7 Ib'^ of CDF Run II data. The 
measured cross sections are well described by NLO pQCD 
predictions including non-perturbative corrections. 
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